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O I We seek to introduce a mathematical method to derive relativistic wave equations for 

two-particle systems. According to this method, if we define stationary wave functions as 
■ special solutions like \I'(ri,r2,i) = ■(/'(ri, r2)e~*^*/'', and properly define the relativistic 

reduced mass ^q, then the new relativistic two-body wave equations can be derived. The 
spin-orbit coupling term is also derived from the relativistic correction terms. On this 
basis, we obtain the relativistic energy levels of pionium and pionic hydrogen atoms, 
using which, we discuss the pair production and annihilation of Tr+ and tt— , and find a 
reasonable explanation for the dark matter. 
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^ ■ 1. Introduction 
o 

Both special relativity and experiments indicate that, the mass of a many-particle 
system in a bound state is less than the sum of the rest mass of every particle forming 
the system, and the difference gives the mass defect of the system, while the product 
of the mass defect and the square of the speed of light gives the binding energy of the 
system. As the binding energy is quantized, the sum of it and the rest mass of every 
particle forming the system is the energy level of the system. For instance, the mass of an 
atomic nucleus is obviously less than the sum of the rest mass of every nucleon forming 
the atomic nucleus. Therefore, in order to express the mass defect explicitly, there is a 
necessity to introduce the concept of system mass, which differs from the sum of the rest 
mass of every particle forming the system. By introducing the concept of the system mass 
and applying proper mathematical skills, the relativistic wave equations for two-particle 
systems is derived. On this basis, let us properly define the relativistic reduced mass to 
further derive the new relativistic two-body wave equations. 
The main results of this paper are expressed as 
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1. The relativistic two-body wave equations 
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2. In the center-of-momentum frame, it is simphfied as 
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Where mo = moi + mo2, E — mc^, and 111,110,11 respectively denote 
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4. For spin 1/2 particles in the central field, the relativistic two-body wave equations 
is expressed as 
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Where Li is the orbital angular momentum of the first particle, and L2 is that of the 
second one. 

5. In the center-of-momentum frame, it is simplified as 
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Where S = Si + S2, L = Li + L2 be the total spin angular momentum and total orbital 
angular momentum operators of the two-particle system respectively, U{r) = U (ri — r2) 
be the interaction energy between particles. 
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2. Relativistic Two-body Wave Equations and Proba- 
bility Meaning of Their Stationary Solutions 

As we know, for a free particle, its energy and momentum are both constant, therefore, 
it is quite natural to assume that its matter wave is a plane wave. According to the de 
Broglie relation 

hk = p, E = hu), 
we have the wave function of a free particle: 

^p(r,t) = v4exp(-z(Et - p ■ r)/^). (1) 

Where E is the total energy of the particle containing the intrinsic energy mgc^, and p is 
the momentum of the particle. 

The quantization method of quantum mechanics is assuming E and p are correspond- 
ingly equivalent to the following two differential operators: 

d 

E -> ih—, p -> -ihV. 
at 

This relation is obviously tenable for the wave function of a free particle, while arbitrary 
wave function is equal to the linear superposition of the plane waves of free particles with 
all possible momentum, namely 




^'(r,t) = / / / c{p,t)'^p{r,t)dp-,dpydp^, 

(2) 
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Where \l'p(r,t) is expressed by ([T]), A = {2txK)^^^'^ . Therefore, this quantization method 
is also tenable for arbitrary wave functions. Let K and K' be two inertial systems, when 
K is transformed to K' , the difference between {E't' — p' ■ v')/h and {Et — p ■ r)/^ is at 
most an integral multiple of 27r, and {Et — p ■ r)//i is a relativistic invariant because of 
the phase periodicity. Then we have 

Aexp[-i(Et - p ■ Y)/h] = Aexp[-i{E't' - p' ■ r')/h] or %(r,i(;) = ^p/(r',t')- 

Thus there is no contradiction between ([2]) and relativity, that is, this quantization method 
itself is relativistic. Clearly, (j2]) are the Fourier transform and its inversion, which can 
be extended to many-particle systems. For two-particle systems, let ri = {xi,yi,zi) 
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and r2 = (0:2, 2/2, -^2) be position vectors of two particles in the laboratory reference 
frame respectively, and corresponding momentum vectors be pi = {Pxi,Pyi,Pzi) and 
P2 = {p 2:2 5 5 ) • Thus related physical quantities in ([1]) and (|2]) are extended to 
r = (ri,r2), p = (pi,P2), A = {27ih)~^, dp^ = dp^^dp^^, dpy = dpy^dpy^, dp^ = dp^-^dp^^, 
dx = dxidx2, dy = dyidy2, dz = dzidz2- Therefore this quantization method can also 
be applied to many-particle systems. Then why the Schrodinger equation derived from 
this method is non-relativistic? That is because the relation between E and p is non- 
relativistic, which is due to the fact that the kinetic energy is non-relativistic. Therefore, 
if we want to establish the relativistic wave equations, we need to introduce the relativistic 
kinetic expression. 

Assuming that any particle with the rest mass mo, no matter how high the speed is, 
no matter it is in a potential field or in free space, and no matter how it interacts with 
other particles, its kinetic energy is: 



This assumption has already been verified (See larXiv: 1008.4224 1 . On this basis, we can 
establish the relation between the system energy E and the momentum p using proper 
mathematical skills, thus obtain the relativistic Hamiltonian. Therefore, we introduce 
the mathematical method in Reference [I]- [9]. Power functions and exponential functions 
play special roles in this method, which are called the base functions as we can establish 
mapping relations between them and arbitrary functions in a certain range. For instance, 
in quantum mechanics, ([2]) determines the mapping relations between wave functions 
of free particles and arbitrary wave functions. Similarly to Reference [1] or abstract 
operators (See a rXiv:1008.3808) . we can introduce the base function and relevant concepts 
in quantum mechanics: 

Definition 1 The right-hand side of ([I]) is defined as the base function of quantum 
mechanics, where E and p are called the characters of base functions, while E and p 
are not only suitable for free particles, but also suitable for any system, and the relation 
between E and p is called the characteristic equation of wave equations. Different system 
has different characteristic equations. 

According to differential laws, we have 



Definition 2 Let mo = moi +mo2 + ■ ■ ■ + moAr be the total rest mass of an A^-particle 
system, E' be the sum of the kinetic energy and potential energy of all the particles, 
then the actual mass of the system, which is called the system mass, is defined as 




(3) 



(4) 



m = mo H — ^-E'. 



(5) 
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Definition 3 If tlie system is in a bound state {E' < 0), tlien tfie absolute value of 
E' is 

I 771/ 1 2 2 A 2 

\E \ = rriQC — mc = L\mc , 

which is called the binding energy of the system, where Am = mo — m is the mass defect 
of the system. 

Definition 4 The total energy of the system E is defined as the sum of the rest 
energy, kinetic energy and potential energy of all the particles forming the system, namely 
E = moc^ + E' 

According to Definition 2 and 4, the total energy of the system is equal to the product 
of the system mass and the square of the speed of light, namely E = mc^, thus the system 
mass is uniquely determined by the energy levels of the system. 

Definition 5 In relativistic quantum mechanics, the stationary wave function for 
two-particle system is defined as the following special solution: 

^(n, r2, t) = ^Pin, ra) exp{-iEt/h). (6) 

Where E is the total energy of the two-particle system. 
Applying (|2]) to Definition 5, we have 

00 
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c(Pi,P2) = (2^/^)3 / f I ^(ri,r2)exp |^-^(pi -ri -Fp2 •r2) ) dxdydz. 
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Where ri = (xi,yi,Zi) and r2 = {x2, 1/2,^2) are position vectors of two particles in the 
laboratory reference frame respectively, and dx = dxidx2, dy = dyidy2, dz = dzidz2- Cor- 
responding momentum vectors are respectively pi = {pa:^,py^,pzj and p2 = {p 
and dp^ = dp^^dp^^, dpy = dpy^dpy^, dp^ = dp^^dp^^. 

In relativistic quantum mechanics, due to the relativistic effect that mass varies with 
speed, the center of mass system is no longer a proper description framework, instead, we 
use the center-of-momentum frame, which is a coordinate system that the total momentum 
equals zero. If fi, V2 respectively denote the speed of particles in the two-particle system, 
then their momentum respectively are 

moiVi mo2V2 
pi = mivi = —====, p2 = m2V2 



and pi = — P2- If V denotes the relative speed between two particles, then we can 
properly define the relativistic reduced mass ^ to make the relative momentum p = /iv 
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satisfy |pi| = |p2| = |p|, namely 

pI, + pI + pI = pI, + pi + pI =pI+pI+ pI (8) 

In other words, the reduced mass can be determined using (jS]) and the relativistic 
velocity addition formula. As it is related to speed, in order to distinguish it from another 
type of reduced mass, we call this one the speed-type reduced mass. For instance, if two 
particles of a two-particle system are restricted to movement along the same line, then its 
speed-type reduced mass is defined as 

H = ■ 1 + ^1. 9 



Where 



moi mo2 

mi = , ^ ==, 7722 



Therefore, using the center-of- momentum frame in (JTj), p2 = — Pi, |pi| = |p|. Substi- 
tuting them into ([7j), we have 
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Where |p| is the relative momentum of the two-particle system. Relative coordinate is 
denoted by r = ri — r2, then the result can be expressed by a more symmetric form 

00 
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Where \E'(ri,r2,t) = ip{ri,r2) exp{—iEt/h) is the relativistic stationary wave functions 
for the two-particle system. Therefore, p = /iv, which is the relative momentum of 
the two particle system in the center-of-momentum frame, is definitely equivalent to the 
differential operator —ihV with respect to the relative coordinate r = ri — r2, still denoted 
by p — )■ —ihV. This quantization rule is tenable for arbitrary stationary wave functions 
for two-particle systems. 
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Considering an isolated two-particle system, if the interaction energy between two 
particles is denoted by f/(ri,r2), then according to Definition 2 and (|3]), we have 

7-1/ TT / 2 2 I 2 4\l/2 2 I / 2 2 I 2 4\l/2 2 ^loN 

Where moi, are the rest mass of two particles respectively, and corresponding mo- 
mentum are pi = |pi|, = |P2|- This is the characteristic equation of relativistic wave 
equations for the two-particle system, thus we obtain. 



This is the spin- less Salpeter equation (See [6], [7]), it is an important relativistic two- 
body wave equation. In order to make it easier to solve the corresponding relativistic wave 
equation, the characteristic equation (fT2l) should be transformed to remove the fractional 
power, then we have 

/771/ TT , 2, 2\2 22, 24,22, 24 

[E ~ U + ruoic + mo2C ) = c mg^c + c P2 + niQ^c 



Expanding the left-hand side of (1141) and applying ([5]), we have 

(mo + m)E' - 2mU + U^/c^ 
= pI+pI + 2{pI + ml,c^y/\pl + ml^c^y/^ - 2moimo2cl (15) 

Further, removing the radical sign, we have 

[(mo + m)E' + 2moimo2C^ - pI ~ pI - "^rnll + U'^/c^f 
= A{pl + ml,c^){pl + ml2^). (16) 

Definition 6 In relativistic quantum mechanics, a type of relativistic reduced mass 
/io of two-particle systems is defined as 

2moimo2 I , 

/iQ = ■ , fx = Ho + -i^E . (17) 

mo + m & 

Where mo = moi -l- mo2, m is the system mass of the two-particle system, E' is the sum 
of kinetic energy and potential energy of the two particles, [i is called the system mass 
corresponding to yUo- Unless otherwise stated, the reduced mass referred in our paper from 
now on is defined in this way, which should not be confused with the speed-type reduced 
mass mentioned previously. 
According to (ITTI) . we have 

(mo + rn)E' + 2moimo2c2 = (mo + m)/ic^ 
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((mo + m)E' + 2moimo2c2)' = (mo + mf{i,, + i^^E' + Aml,ml^c\ 
Thus (fT6|) can be expressed as 

/l/22| 2 22, 2 2 2\ 

= (mo + m)^(/io + /i)c^E' - 2(mo + m)^c^{j)l + pi + 2mU - U'^/c^) 
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According to f|T8l) . the relativistic Hamiltonian of two-particle systems can be ex- 
pressed as 
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Therefore, taking (fT9|) as the characteristic equation, multiplying both sides of the equa- 
tion by the base function \E'p(r, t), and by using (jl]), we have 
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Where f/(ri,r2) denotes the potential energy of the interaction between two particles, 
V^, V2 are Laplace operators respectively corresponding to ri,r2. 

According to ([2]), in the operator equation which is tenable for the base function 
\E'p(r,t), as long as each operator in the operator equation is a linear operator and each 
linear operator does not explicitly contain the characters E and p of \E'p(r,t), then this 
operator equation is also tenable for an arbitrary wave function \E'(r, t). Whereas, consid- 
ering that the system mass m is equivalent to the character E, this operator equation is 
not tenable for arbitrary wave functions, but tenable for an stationary wave function like 
([6]). Thus we have: 

An isolated two-particle system, the total spin angular momentum of which is zero, 
is described by the stationary wave function \l'(ri, r2, t) or '?/^(ri,r2), any stationary wave 
function 

^'(ri, r2, t) = ^(n, r2) exp{-iEt/h) 
satisfies the following relativistic wave equation 
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Where mo = moi + mo2, E' = E — rrioc^, E = mc^. m, fio, fi respectively denote 

1 , 2moimo2 1 „/ 

m = mo + —E , fio = ■ , fi = no + —E . 

mo + m 

Clearly, for non-relativistic limits, we have 

moimo2 moimo2 



/i /iQ 



mo moi + mo2 



In other words, the relativistic wave function ilj{ri,r2) for two-particle systems is 
determined by the following relativistic wave equation and natural boundary conditions: 
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For bound states, the total energy E of the system is quantized, which is called the 
system energy level. The system mass m = E/c^ is uniquely determined by the system 
energy level E. Clearly, for non-relativistic limits, this equation turns out to be the 
Schrodinger equation of two-particle systems. If the system is in the external field, then 
the system potential energy U{yi, Y2) includes both the potential energy of the system in 
the external field and the interaction energy between particles. 

Using the center-of-momentum frame, then according to ([8]) we have p\ = P2 = P^i 
where p is the relative momentum. Considering 
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then in the center-of-momentum frame, (fT9l) is simplified as 
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Taking (123!) as the characteristic equation, similarly we have: Considering an isolated 
two-particle system in the center-of-momentum frame, if the total spin angular momentum 




11 



of the system is zero, then the stationary wave function 

is determined by the following relativistic wave equation and natural boundary conditions: 
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These are also the expressions of relativistic wave equations (|20l) and (!2Ti) in the center- 
of-momentum frame respectively, where is the Laplace operator corresponding to the 
relative coordinate r = ri — r2. 

Relativistic wave equations and fl2^ can be further expressed as 
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If ^(r,t) is the solution of ([26]), then for M;(r,t) = ^*(r, t)^(r, t) we have 



9ty 



2m 



ih 



+ 



mo + m / iJio + II 
1_ ^ 

ih (mo + m)2(/io + /i)c^ 
1 4/1^ 
ih (mo + my{fiQ + ji)c^ 

2 



2mU - ^ ) ($*V^^ - W^^*' 



2m 



mo + mj /iQ + /i 

42^ 



V ■ (**V^ - ^V^ 



Vf 2mf/ 



(mo + m)2(/io + /i)c2 \ 

Aih ( IP^ 

(mo + m)^(/io + /ijc^ V 



2m 



mQ + mJ /io + 
(mo + m)2(/io + /i)c^ 



V ■ (\I/*V\[^ - \I>V^* 



V- 



2mf/ - — ) (^*V^ - ^V^*) 



(27) 



In other words, according to (!26l) . w(r,t) = \l/*(r, t)\l/(r, t) is defined as the probability 
density, then the corresponding probability current density vector J is expressed as 



2m 



ih 



Aih 



(mo + m)2(/io + yu)c2 V 
Aih 



2mU 



(mo + my{j2o + jj,) 



im^ - 2m— + — j ($V** - ^*V^). 
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Therefore, similar to non-relativistic quantum mechanics, for a two-particle system 
that the total spin angular momentum of the system is zero, we have 



dw 
'di 



Where J is expressed as 
J 



+ V- J = 0. 



u 



(2^ 



(29) 



(mo + m)2(^o + Ai) V 

On the right-hand side of fl27|) . if combining the first and the fifth terms, the stationary 
relativistic wave equation for two-particle systems in the center-of-momentum frame is 
expressed as 



E'lp 



4r 



u 



+ 
+ 
+ 



(mo -I- m)2(/io + /i) 
2/i 

(yUo + Ai)(mo + m) 
{mo + m)2(/io + /i)c2 



m 



2mU - — 1 ^ 



( 



2mU - 



(mo + mY{^Q + /i)c"= \ / 
1 / ^ f/^Y 

(mo + m)2(/io + /i)c2 \ c2 y 



(30) 



3. Relativistic Two-body Wave Equations and Spin 
Angular Momentum 

As we know, the spin angular momentum not only has the general property of angular 
momentum S x S = i^S, but also has its own particularity. For instance, for electrons 
and protons, the projection of the spin angular momentum S in any direction only takes 
two values ±/i/2. 

For the convenience of studying this type of angular momentum, a type of dimension- 
less vector a is introduced, determined by 



(T X (J = 2icr, = Gy = a1 = 1. 



(31) 



Using this type of vector a, the spin angular momentum is expressed as S = {h/2)a. 
According to (|3T1) . if A, B are two arbitrary vectors commuted with a, we have 



((T ■ A)((T ■ B) = A ■ B + i(T • (A X B). 



(32) 
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If the two particles composing the system both have spin 1/2 angular momentum, 
denoted by Si and S2: 

J5l = -(Ti, ^2 = -(72, 

then the relativistic wave function describing the system should be a spinor. Let 0, ip be 
arbitrary spinors, then we can easily prove: 

(a ■ p)(0t^) = 0t(^ . p)^ + . p0t)^. (33) 
Repeatedly using ( 15^ . considering (a ■ p)((T ■ p) = p ■ p = p^, we have 

p2(0t^) = ^ 2(a . p(/)^)((T ■ p^) + (pV^)V'- (34) 

Therefore, under the condition that wave functions are spinors, the relativistic Hamilto- 
nian of two-particle systems f ll9p is rewritten as 



E' 



2(moi/i + mo2/io) pj ^ 2(mo2/i + moi/ip) 
(mo + m) (/io + /i) moi (mo + m) (/io + /i) mo2 



2/i / 2m ^ f/2 
2 

(mo + m)2(/io + 



(Ji • Pi I 2mU - 
CT2 ■ P2 I 2mf/ 



{nio + mY{^Q + /i)c2 

? 2mU - — ^ (pi + pI) 

(mo + m)2(/io + /i)c2 \ J ^ ^ 

(p2 _^ f _ 

(mo + m)2(/io + /i)c2 [ ^ ^ V y 

1 / 2m _ \ ^ 

(/io + /i)c^ \mo + m {mo + 171)0^ J 

(mo + m)2(/io + /i)c2''^^ 



(f^i ■ Pi) 
(0-2 • P2) 



(35) 



Thus the stationary relativistic wave equation for two-particle systems can be expressed 
as 
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2(moiAi + mo2/io) Z?-^ ^2,j, 2(mo2/i + mpi/io) /i^ ^^^^^ 



(mo + m)(yUo + Ai) 'tioi 



(mo + m)(/io + A*) 



+ 



2/i 



2m 



-f/ 



f/2 



+ 



+ 



+ 



+ 



/io + yU \mo + m 

(mo + m)2(/io + /i)c2 
2/1^ 

(mo + m)2(/io + /i)c2 
(mo + m)2(/io + /i)c2 
(mo + m)2(/io + yu)c2 



(mo + m)c2 



(72 ■ V2 

2m[/ - 

(V? + v: 







|^2mf/ - 




c2 ) 






(^2mU - 




C2 J_ 



(ai ■ Vi)^ 

((T2 ■ V2)^ 



2m 



-U 



2mU 
2 



LP 

C2 



(/io + /i)c2 ymo + m (mo + m)c2 

^ ^ . , (V^ - V^)2vl/ + moc^v^. 

(mo + m)^(/io + /i)c^ 



(36) 



The two particles of the system are spin 1/2 particles, where U{ri, r2) denotes the potential 
energy of the system in the external field and the interaction energy between particles. 

In the central field, the potential energy U{r) is only related to the distance between 
particles, but not related to the direction of r. According to ( 132|) . for the fourth and the 
fifth terms on the right-hand side of (136|) . we have 



2h' 



(mo + m)2(/iQ + /i)c2 
2h' 

(mo + m)2(/io + /i)c2 



f/2 

a ■ V 2mU ^ 



(a- V^) 



V ( 2mU - ^ 1 ■ 



+ 



2r 



(mo + m)2(/io + /i)c2 
(mo + m)2(/io + /i)c2 



(mo + m)2(/io + /i)c2 



m 



m — 



V I 2mU ^ 1 X 

U\ dUd^ 



S • L*. 



c^y dr dr 
U\ IdU 



J r dr 



Thus in the central field we have 

For spin 1/2 particles, the stationary wave equation for two-particle systems can be 
express as 
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(mo + m)(/io + Ai) moi 
2/i / 2m 



f/2 



(mo + m)(//o + A*) mo2 



+ 
+ 
+ 



4/^2 (m - U/c^) / dt/ 9^ dU (9^ \ 
(mo + mY{iJio + /x)c^ 9ri dr2 dr2 J 



(mo + m)2(/io + /i)c2 

(mo + m)2(/xo + /x)c2 
1 / 2m 



2mU 



(V? + V^)* 



(V? + V^) (^2mU - — 
U--, r^l * 



{/jLo + IjI')c'^ \mo + m {mo + m)c 

- U/<^') • L,* + 1^S2 • 

(mo + my{iiQ + ii)c^ \ri dri r2 dr2 



(mo + m)2(/xo + //)c2 



(Vt - V^)'* + moc^*. 



E'lp = 



2(moi/x + mo2Ato) ^2^, _ 2(mo2At + moiAto) ^2^, 
(mo + m) (/xo + A«) moi ^ (mo + m) {/lo + jj) mo2 ^ 



+ 



+ 



2/i / 2m [/^ 
IJ,o + H \mo + m (mo + m)c^ / 



m 



fdUdi) dU dip 
(mo + m)^(/io + A*)c^ \(iri 9ri (ir2 (?r2 / 

, ^ ^ , , f 2mC/ - ^) (V? + V^)^ 

(mo + m)^(/xo + Aijc^ V <^ / 



+ 



+ 



(mo + m)2(^o + 



iVl + Vl)l2mU-^ 



{Hq + /x)c2 \mo + m ^ (mo + m)c2 / ^ 

8(m-C//c2) /IrfC/^ ^ , IdU^ ^ \ 

(mo + my [1x0 + ix)c^ \ri dri r2 dr2 J 

t-. 1 



(37) 



(38) 



(mo + m)2(//o + //)c2^^ "^^J V^- 

Where Li is the orbital angular momentum of the first particle, and L2 is that of the 
second one. ^ is the stationary relativistic wave function 

^(ri, r2, S12, S2z, t) = V^(ri, r2, si^, S22) exp{-iEt/h). 
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For an isolated two-particle system, in the center-of- momentum frame, p2 = — Pi, |pi | = 
|p|, and IpI is the relative momentum of the two-particle system. Based on the corre- 
sponding relation between momentum operators and gradient operators, according to 
P2 = — pi = — p we have V2 = — Vi = —V. Where Vi, V2 and V are gradient operators 
corresponding to the coordinates ri, T2 and r = ri — r2. 

Therefore, if using the center-of-momentum frame in the wave equation (!37|) . suppos- 
ing an isolated two-particle system, then Using V2 = — Vi = —V and p2 = — pi = — p, 
in the central field we have 

dU d<i/ „^ „^ dUd<I/ 

— — = (Vit/ ■ Vi^ = (Vf/) ■ = — — . 
ari OTi ar or 

dU „ ^ , . dUd^ 

-^t:- = V2f/ ■ V2^ = i-VU) ■ -V^ = — — . 
ar2 or 2 dr or 

1 dU^ 1 dU IdU 
— — Li = — :^ri X pi = (Vit/) X pi = [VU) X p = -— r x p. 
ri dri ri dri r dr 

1 dU^ 1 dU ,^ , , , IdU 

— — L2 = — r2 X p2 = {V2U) X p2 = (-Vf/) X (-p) = -— r x p. 
r2 dr2 r2 dr2 r dr 

Considering (!22|) . the wave equation (137|) can be expressed as 



2fe2 



2m 



U 



f/2 



+ 



(/io + /i)(mo + m)2 /io + V"^o + ("^0 + "^)c^ 

(mo + m)^(/io + Ai)c^ dr dr 



+ 



+ 



4r 



(mo + m)2(/io + /i)c2 



2r 



(mo m,)2(/io + /i)c2 
1 /2m 



2mU - — \ V^^ 



V' ( 2m[/ - ^ 



(Ho + 11)0^ \mo + m (mo + m)c^ 
8(m - U/c^) 1 rft/ 



(mo -|- my{^Q + /i)c^ r dr 



Si + S2) ■ (r X p)^ + moC^^. 



(39) 



The total spin angular momentum of the system is S = Si + S2, and the orbital angular 
momentum L is 



L = Li + L2 = ri X pi + r2 X p2 = (ri - r2) X pi = r X p. 

Therefore, the total orbital angular momentum L of the two-particle system in the center- 
of-momentum frame is equal to the cross product of the relative coordinate r and the 
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relative momentum p. Combining the first and the fourth terms on the right-liand side 
of fl39|) . we have 

Let S, L be the total spin angular momentum and total orbital angular momentum 
operators of the two-particle system respectively, U{r) be the interaction energy between 
particles, then in the central field, the stationary wave function for the two-particle system 
in the center-of-momentum frame 

*(r, s^, t) = ipir, s^) exp{-iEt/h) 

is determined by the following relativistic wave function and natural boundary conditions, 
namely 



ih 



(9* 
'dt 



4n' 



(/io + /i)(mo + m)2 
2/i 



m 



U 



+ 
+ 



2mU ^ 1 ^ 



(Ho + l^){mo + m 
(mo + m)^(/io + yu)c^ dr dr 



2r 



(mo + m 


)2(/io + /^)c2 
1 


{niQ + m 




8(m 





V2 ( 2mt/ - ^ 



2mU - 



U 



2\ 2 



(mo + m)2(/ig + /i)c^ r dr 



— — S ■ + moc^^'. 



(40) 



4h' 



(/io + /i)(mo + m)2 
2/i 



U 



m — 



+ 



(/io + /u)(mo + m) 



2mU -]^p 



+ 



+ 



8n\m-U/c^) dUdtjj 
(mo + m)^(/io + /i)c^ (ir 9r 
2/1^ 



( 



(mo + my{fio + /i)c2 
1 

(mo + m)2(/io + /i)c2 \^ 

8(m - U/c^) 1 rff/ 
(mo + m)2(/iQ + /i)c^ r rfr 



f/2 

2mt/- — 



2mf/ - 



S ■ Lt/-. 



(41) 
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According to the wave equation (l30l) . in the central field, the stationary relativistic 
wave equation for the two-particle system that the total spin angular momentum is zero 
in the center-of-momentum frame is expressed as 



Ah' 



+ 



+ 



+ 



(mo + m)2(^o + /u) 
2/i 

(/io + /i)(mo + m) 
(mo + m)2(/io + /i)c2 

(mo + m)2(/io + fJ,)c^ 
1 

(mo + my{fio + /i)c2 



m \ w %p 



2mU 



( 2mU - — 



U\ dUdtp 
m 



J dr dr 
2mU t) i^- 



(42) 



Compared with ( H2|) . (HTj) has one more term, corresponding to the coupling energy be- 
tween the total spin angular momentum and the total orbital angular momentum. 

In (I35l)-(l42l). E' = E — niQC^, mo = moi + mo2, E = mc^, m, fiQ, /i respectively denote 

1 , 2moimo2 1 „/ 

m = mo + —E, /io = ■ , fi = fio + ^E . 



If applying the relativistic wave equation ( 1411) for two-particle systems to hydrogen-like 
atoms, then U/ ^ m, (141 p can be approximately expressed as 



^, , 4m^ „n , ArnaU , 4m^ f/^ 

+ , ^fft^^^^_i^|^S.L^. + 5!(V^wV (43) 



(mo + m)^(/io + lJ,)c^ \ dr dr r dr 2 

Applying this wave equation, we expect to determine the hyperfine structure of the energy 
levels of hydrogen-like atoms and the Lamb shift with high precision. 



4. Relativistic Energy Levels for Two- Particle Systems 
with Zero Total Spin Angular Momentum 

A hydrogen-like atom that the total spin angular momentum is zero, and the pionium 
composed by 7t~ and tt"*", are both two-particle systems. This type of potential energy of 
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interaction between particles is U = —Zej./r, = e(47r£:o) considering 

r 

the wave equation P2l) can be further expressed as 



Ah' 



m + 



c^r 



2\ 2 



(mo + m)2(/io + yu) 
(/io + /i)(mo + m) V 



+ 



^2 1^ 2 



8r 



m + 



1 /2mZe? Z^elA' 



(mo + m)2(yUo + /u)c^ \ r c^r^ 



2„4 



(mo + m)^(/io + yu)c^ c^r"^ 



{mo + m)2(/io + yu)c^ 



(5(r)'?/^. 



(44) 



Now let us solve the wave equation ( l44l) under the condition of r > 0, when 6{r) 
Using the spherical polar coordinates, the Laplace operator is expressed as 



dr 



dr 



sine 89 V'' 89 



8^ 



sin 9 8ip'^ 



Referring to the solving procedures of the Schrodinger equation in Ref . [5] , supposing 
'ip{r,9,{p) = R{r)Y{9,ip), substituting it into (jH]) and considering S{r) = 0, we have 



'^2 



(mo + rriY^HQ + /i)_EV 



1 d 



,dR 



+ — 

+ 



2Ze' 



1 dR 



Rdr \ dr J [m + Zel/ {c^r)]c^ R dr 
fj,{mQ + m)r^ ^ 2mZe^ ^ Z^ej 

^2 ^ 2 



2r[m + Ze2/(cV)]2 V r 



+ 

+ 



4^^c2[m + Ze2/(c2r)]2 V r 
1 Z^e! 



+ 



^2^2 



[m + Ze^/(c2r)]2 cV^ 



1 

T 



1 9 / . 8Y\ 
sm9—-— 



sin 6' 



1 8^Y 



89 J sin^ 9 8ip'^ 



A. 
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+ 



mc^r J dr 
IxirriQ + m)Ze 



11 (^2dR\ _ , ^ 



1 + 



mh^r 
Zei 



2 



Zel \ 2Zel dR 



R 



mc^r J mc^r"^ dr 
{mo + m)^(yUo + fJ')E' 



Am?h 



R 



mc^r j 



A 



i? = 0. 



(45) 



1 d ( . dY_ 



1 d^Y 
sin^ 9 d(p'^ 



XY = 0. 



(46) 



According to fH6|) . denoting A = Z(/ + 1), / = 0, 1,2, . . ., obviously the solution of the 
equation is the spherical harmonics Yim{6, 

Now let us solve the radial equation fH5|) . discussing the situation of the bound state 
{E' < 0). Let 



«' = ^^^^[(/io + /i)|i?1]^^ P = «V, (47) 



(3 



mh 



1/2 



a'm/i^ [(/io + /^)|^'L 

using the variable substitution p = a'r, then f HS]) can be expressed as: 



(4J 



(49) 



Where a denotes the fine structure constant, do, which is a small parameter, denotes 



do = 2Z^a'D, D 
Let -R(p) = u{p)/p, considering 



p(mo + m) 
2m2 



1 d f2dR\ _ 1 
p^ dp \ dp J p dp^ ' 



(50) 



then fl49l) can be expressed as: 



(in 



2 ^2. 



^ Pp) dp^ \ ^ (3pJ Pp^dp^ Pp^""^ p 



1 ^ A ^ ^0 



/3p 



2/3p 
/(/ + !) 



0. 



(51) 
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Firstly, let us study the asymptotic behavior of this equation, when p — )■ oo, the 
equation can be transformed into the following form: 

(fu 1 

— - -u = 0, u{p) = exp(±p/2). 

As exp(p/2) is in conflict with the finite conditions of wave functions, we substitute 
m(p) = exp(— p/2)/(p) into the equation, then we have the equation satisfied by /(p): 



dn\'^ <ff ( dn\ / do 2d() \ d f 



l3p J dp^ \ /3p J \ /3p Pp"^ J dp 



Thus solving for the radial wave function -R(p) comes down to solving for /(p), namely 

, , 1 / 1 \ , , 2Z 2m 

i? p = -exp --P / p , p = — r, ao = ■ ^. 53 

p \ 2 J pao rriQ + m pej 



According to (llSil we have 

p2 _ P 



{po + p)\E'\' 

Substituting p = pq — \E'\/c^ into the equation above, we have 

(Z^a^ + - 2poc'(^'«' + I3')\E'\ + plc^Z'^a' = 0. (54) 

Solving fl5^ . we obtain two roots of \E' 



p'l- 



l^'l = poc^^po& 1 



-1/2 



Thus we obtain the system mass p corresponding to the reduced mass po^ 

p = ±po M + ^1 . (55) 

According to Definition 6, po, respectively denote 

2moimo2 1 , 

A^o = : , fJ' = fJ'0 ^1^1- 

mo + m 
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Considering m = mo — \E'\/c^, mo = moi + mo2, \E'\ is derived by fl55|) . As the total 
energy of the system E = tjiqC^ — \E'\ = mc^, we can further obtain the total energy E 
and the system mass m. According to ( !55l) . fi has positive and negative values. When /i 
takes on positive values, E is expressed as 



E = ± 



niQ^ + 2moimo2 1 + 



2^,2 \ -1/2 



m, 



02 



1/2 



When /i takes on negative values, E is expressed as 



E = ± 



m 



01 



- 2moimo2 1 



^1/2 



m; 



02 



1/2 



Therefore we have two positive and two negative energy solutions. Negative energy solu- 
tions are related to the ubiquity of antimatter, which will not be discussed in this paper. 
Taking on positive energy solutions, we have 

1/2 

c^. (56) 



E 



moi + 



2moimo2 



E 



m, 



2moimo2 



+ "^02 



01 



+ m^2 



1/2 



(57) 



Corresponding to the two positive energy solutions, the system mass m has two ex- 
pressions: 

2moimo2 



m 



m 



moi + 



mn 



2moimo2 



+ m^2 



1/2 



1 1/2 



+ m. 



02 



(58) 



(59) 



Clearly, when (!55l) takes on positive values, the system mass m is expressed by (!58|) . But 
when (l55l) takes on negative values, m is expressed by (!59|) . 

In (1561) and (J57l) . the quantization of E is mirrored by the fact that (3 is related to both 
the principal quantum number n and the angular quantum number I, i.e. P = P{n,l). 
Solving the equation fl52l) we obtain the expression of f3{n, I). Therefore fl56l) and (l57j) are 
the general expressions of the relativistic energy levels for two-particle systems. 

It seems difficult to accurately solve (l52l) . Let us solve this equation for approximate 
solutions to obtain the approximate expression of f3{n,l). First, (l52i) is expressed by the 
standard form of second-order ordinary differential equations 



f" + pip)f' + qip)f = 0. 



24 



Where p{p), q{p) respectively denote 



P{p) 



1 + 



do 

Pp 



-1 



T _|_ ^ _|_ ^'^0 

Tp J? 



<l{p) 



dn 



+ 11 + ^1 U^a. 



do \ do 1 



+ 1 



/3p 



-2 



{Z'a' + I) + 12] 



dl 1 



4/32 p^' 

Considering do is very small, we have p{p) ~ —1, and q{p) is approximately expressed as 

c^o 1 1 I „o o 3do do\ 1 



( !52l) is approximately expressed as 



(ip2 



^ + 



2PJ P 



2^2 



/ = 0. (60) 



Thus p = is a regular singular point of the equation (160|1 . Suppose the series solution 
of this equation can be expressed as 



(61) 



u=0 



In order to guarantee the finiteness oi R = u/p ai p = 0, s should be no less than 1. 
By substituting (l6T]l into (160|) . as the coefficient of p'^'^'^~^ is equal to zero, we have the 
relation satisfied by b,^: 



bu 



+1 



s + v-[P + do/{2P)] 



{s + p){s + p + + + Z^a^ - 3do/2 + do/ 13 



b.. 



(62) 



If the series are infinite series, then when i/ — )■ 00 we have b^j^i/b^ — )■ 1/z/. Therefore, 
when p — )■ 00, the behaviour of the series is the same as that of e'', thus /(p) in ( l53l) tends 
to infinity when p — )■ 00, which is in confiict with the finite conditions of wave functions. 
Therefore, the series should only have finite terms. Let 6n,.P*^"' be the highest-order 
term, then bn^+i = 0. By substituting 1^ = 71^ into ( 162|) we have (3 + do/ (2/3) = n,, + s. On 
the other hand, the series starts from 1^ = 0, therefore, 6_i = 0. Substituting u = —1 into 
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fl62|) . considering 60 7^ 0, we have s{s — 1) = /(/ + 1) — Z'^a'^ + Mo/2 — do/13. Denoting 
n = rir + I + 1, then the following set of equations can be solved for s and (3: 



s{s - 1) = /(/ + 1) - ^2^2 + ^do/2 - do/P 
/3 + do/{2P) =nr + s 

n = Ur + I + 1 



(63) 



We derive s = 1/2 ± ^{l + 1/2)2 _ ^2^2 + 2>do/2 - do/ (5. Considering s should not be 
less than 1, s = 1/2 + a/(/ + 1/2)2 _ z'^a'^ + 3do/2 — do/ (3. Thus we obtain a specific 
expression of (3{n, I): 



do_ 

2/3 



+ \ 1 + 



-Z^a^ + ^-^-^ = n-a, (64) 



Where a, = / + 1/2 + do/i2l3) - ^{l + 1/2)2 _ ^2^2 + 3^^/2 - do/ (3. 

Therefore, the relativistic energy levels for two-particle systems (l56l)-( j57l) . the system 
mass (I58l)-(l59l) and (l55l) can be respectively expressed as 



Er, 



niQ^ + 2moimo2 1 + 



Z2«2 



mgi - 2moimo2 ( 1 + 



m 



m 



niQ^ + 2moimo2 1 + 



01 - 2moimo2 ( 1 + 
Z2« 



= ±/io 1 + 



{n - aiY 

{n - aiY 

{n - aiY 

{n - aif 
2^,2 \ -1/2 



1 1/2 



-1/2 



+ "^02 



-1/2 



-1/2 



1/2 



1/2 



1/2 



+ "^02 



n - o-/ 



2moimo2 
mo + m 



1 



dn 



2 2(n - aA 



1 + 



do = 2Z'^a^D, D 



- Z^a^ + 



yu(mo + m) 
2m2 



3do 



do 



n~ ai 



(65) 

(66) 

(67) 

(68) 

(69) 

(70) 
(71) 



Considering a two-particle system, which is a hydrogen-like atom composed by a spin- 
zero nucleus (like the deuteron) and a tt", called a pionic hydrogen atom, moi, mo2 are 
the rest mass of vr" and nucleus (moi ^ ^^02), then ( 165|) can be expanded as the following 
fast convergent infinite series 
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Clearly, we obtain the normal energy levels for two-particle systems. Using fl65|) . we can 
calculate the energy spectrum of pionic hydrogen atoms more accurately. The energy 
levels expressed by fl66l) are called the abnormal energy levels. Unlike the normal energy 
levels, the abnormal energy levels decrease with increasing the principal quantum number 
n. 

According to fl65l) - (f?T!) . we need to use iterative methods for calculation. As is very 
small, taking = in the expression of o"/, we can obtain the zeroth order approximation 
of (7/ to calculate that of m and /i, then substitute them into do to obtain its zeroth 
order approximation. Substituting the zeroth-order approximation of d^ and ai into ai 
to calculate its first-order approximation, repeating this process, using the first-order 
approximation of ai to calculate that of m and /x, substitute them into rfg to obtain 
its first-order approximation. Substituting the first-order approximation of d^ and ai 
into ai to calculate its second-order approximation, and repeating this process we can 
calculate the nth-order approximation of o";, further, we can obtain the energy levels En 
and reach the required accuracy. This calculation process can also be realized by computer 
programming. Note that (!69|) means /i has both positive and negative values. When 
calculating normal energy levels, /i > 0, the system mass uses (1671) . When calculating 
abnormal energy levels, fi < 0, m uses (l68i) . Then what is the physical meaning of 
abnormal energy levels? 

A bound state composed of a positive and a negative pion is called the pionium. The 
positronium, pionium, protonium, neutronium, etc., are generally called the particleium. 
For this type of system, we have Z = 1 and moi = Tno2- Therefore, when the pionium is 
at abnormal energy levels, according to ( 166|) we have 



What is the physical meaning of this result? According to En = mc^, we have m = 0, 
which means the disappearance of the particle system and the annihilation of a pair 
of positive and negative pions. In relativistic quantum mechanics, the meaning of the 
vacuum state should not be restricted to a state that the energy is zero. For any bound 
state composed of a particle-antiparticle pair, if it is at abnormal energy levels, it is in the 
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(moi - mo2)c^ = 0. 
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vacuum state. Therefore, the annihilation of a pair of positive and negative pions has two 
phases. The first one composes the pionium, while the second one is its transition from 
normal energy level expressed by (l65l) to abnormal energy levels expressed by ( !66l) . If this 
process produces 7 photons, it means a pair of positive and negative pions annihilates 
into photons. The reverse process is the pair production of positive and negative pions. 
A reasonable extension of this concept is that after the annihilation of any particle- 
antiparticle pair, a small percentage of the energy is generally given to the abnormal 
energy levels of the particleium. Thus this percentage of energy is also quantized, and its 
energy spectrum is given by the abnormal energy levels of the particleium. For instance, 
let rTT-Tr be the rest mass of tt^, then the abnormal energy levels of the pionium can be 
expressed as 



The non-relativistic approximation of this formula can be easily obtained: 



which is completely different from the energy spectrum of normal matter expressed by 
( 165|) . Thus generally, photons emitted by normal matter do not interact with the par- 
ticleium at abnormal energy levels. Assuming there is all sorts of the particleium at 
abnormal energy levels everywhere in the universe, they do not absorb photons emitted 
by celestial bodies because of the huge difference in their energy spectrum. On the other 
hand, as antiparticles have opposite properties from particles, the particleium is equiva- 
lent to a neutral particle with small cross section. When the particleium annihilates and 
transits to abnormal energy levels, this feature of small cross section is reserved. In other 
words, there may be full of energy left behind after the annihilation of particle-antiparticle 
pairs in the universe, each type of energy has its own spectrum and each spectrum is 
determined by abnormal energy levels of the corresponding particleium. Based on the 
relativistic mass-energy relation, the particleium at abnormal energy levels has mass and 
can cause gravitational effects. Apparently, the particleium can transit from higher ab- 
normal energy levels to lower ones with the emission of photons, but this type of photons 
can hardly reach our observation instruments as they will be easily absorbed by other 
particleium at abnormal energy levels nearby, and that is why it seems dark. According 
to observation results of modern astronomy, there exists huge amount of mysterious dark 
energy and dark matter in the universe. If the particleium at abnormal energy levels is 
the carrier of dark energy, and part of them aggregate to form dark matter because of 
gravitational effects, then the abnormal energy levels can establish the theoretical basis 
of the existence of dark energy and dark matter. 




28 



5. Conclusion 



In conclusion, by introducing Definition 1-6, using (jj]) and assuming the relativistic 
kinetic expression is tenable on a wider scale, the relativistic wave equations for two- 
particle systems are derived, and the new relativistic two-body wave equations are ob- 
tained, from which, the spin-orbit coupling term is also derived. By applying this type of 
wave equations to pionium and pionic hydrogen atoms, the general expression and spe- 
cific calculation formulas of relativistic energy levels for two-particle systems are derived. 
Besides, we further find the relativistic abnormal energy levels, thus the pair production 
and annihilation of particles and antiparticles boil down to the transition between normal 
and abnormal energy levels of two-particle systems. Then the relation between abnormal 
energy levels and dark energy is discussed. We believe there are all sorts of particleium 
at abnormal energy levels everywhere in the universe, which are not only dark, but also 
hardly interact with celestial bodies and photons passing through them. Proposing the 
abnormal energy levels paves a new road to reveal the mysteries of dark energy and dark 
matter. 
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